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By considering a simplified but exact model for realizing the ekpyrotic scenario, we clarify various 
assumptions that have been used in the literature. In particular, we discuss the new ekpyrotic 
prescription for passing the perturbations through the singularity which we show to provide a 
spectrum depending on a non physical normalization function. We also show that this prescription 
does not reproduce the exact result for a sharp transition. Then, more generally, we demonstrate 
that, in the only case where a bounce can be obtained in Einstein General Relativity without 
facing singularities and/or violation of the standard energy conditions, the bounce cannot be made 
arbitrarily short. This contrasts with the standard (inflationary) situation where the transition 
between two eras with different values of the equation of state can be considered as instantaneous. 

We then argue that the usually conserved quantities are not constant on a typical bounce time scale. 

Finally, we also examine the case of a test scalar field (or gravitational waves) where similar results 
are obtained. We conclude that the full dynamical equations of the underlying theory should be 
solved in a non singular case before any conclusion can be drawn. 

PACS numbers: 98.80.Hw, 98.80.Cq 


I. INTRODUCTION 

Modern ideas of particle physics, such as superstring, 
M— theory [jl| or quantum gravity ||], cannot in general 
be subject to experimental constraints because of the 
enormous energies (usually of the order of the Planck 
mass) at which they are supposed to become effective. 
According to recent theoretical developments i! , there 
is hope that space possesses more than three large dimen¬ 
sions and that these extra dimensions might turn out to 
be observable in a not too distant future. The majority 
of the theoretical models that have been built so far are 
however based on extremely high energy extensions of 
the standard particle physics model, and thus currently 
need to be tested by the yardstick of cosmology, the latter 
being the only playground at which those theories could 
have acted. 

According to the now standard paradigm that de¬ 
scribes the early universe and that is expected to stem 
from such high energy particle models, a phase of su¬ 
perluminal accelerated expansion known as inflation |^] 
preceded the radiation-dominated epoch. Up to now, no 
model has come as close to being a reasonable challenger 
to solve the standard cosmological puzzles (flatness, ho¬ 
mogeneity and monopole excess). The extra bonus pro¬ 
vided by the inflationary phase is that it leads naturally 
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to a scale-invariant density fluctuation spectrum that 
seems to be in agreement with the observations. 

Inspired by the recent developments of M— theory |j), 
in particular through Re f. |7j l , and invoking brane cos¬ 
mology, recent work [|| ||, |l(| claimed to be able to solve 
all the aforementioned problems as well, including a new 
way of producing primordial cosmological perturbations. 
Although the model, both in its “old” || and “new” || 
versions is plagued with many difficulties JTfij, as a po¬ 
tential alternative to inflation (see also Ref. [|12| in that 
respect), it is worth examining in detail, would it be only 
to re-enforce the confidence we may have in the latter. 


In both the original and most recent versions, the uni¬ 
verse is supposed to consist of a four dimensional (visi¬ 
ble) brane evolving in a higher (in practice 5) dimensional 
bulk. By assuming the brane to be a Bogomolnyi-Prasad- 
Sommerfield (BPS) state |13|, one ensures that the cur¬ 
vature K, vanishes, thus addressing the flatness problem. 
To begin with, another brane, that can be either a light 
bulk brane || , or the other (hidden) boundary brane 0, 
moves freely in the bulk until it collides with the visi¬ 
ble brane. The collision time is interpreted as the hot 
big bang at which point the model is made to match the 
standard cosmological model. 


Apart from the collision time, the theory, which can 
be seen as effectively four dimensional in the long wave¬ 
length limit, relies on the General Relativity (GR) the¬ 
ory together with some extra fields. In this effective 4D 
model, the Universe collapses, experiences a bounce at 
some instant in time, and starts expanding. As far as 
cosmological perturbations are concerned, only GR cal- 
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culations have been discussed up to now. 

The pre-impact phase has been the subject of many 
tentative calculations of the perturbation spectrum 
that would be generated by quantum perturbations of 
the brane [|10| . A general agreement has now been 
reached |l0| |14|, 15, [l(| that the curvature perturbation 


spectrum has spectral index n £ = 3, while that of 
the Bardeen potential P$ ends up with n s = 1, i.e., a 
scale invariant spectrum. On the other hand, the spec¬ 
tra of $ and £ are identical in the post-impact phase, 
and enter the Cosmic Microwave Background Radiation 
(CMBR) multipole moments. It is therefore of utmost 
importance to obtain full knowledge of these spectra not 
only in the pre-impact phase, but also after the bounce 
has occurred, i.e., at times that are observable now. In 
other words, the fate of <f> and £ through the bounce is 
the main issue before any conclusion regarding the model 
can be drawn. 


Only a few definite statements can be done about the 
bounce epoch. The first, which was advocated by many 
authors, is that GR does hold during it, or, stated dif¬ 
ferently, that it lasts sufficiently little that corrections to 
GR can be regarded as negligible. Lacking the actual 
theory, this is the only statement that can be endowed 
with a predictive power. To begin with, it implies that 
there was no singularity, and, if the null energy condi¬ 
tion is to be satisfied, that space is positively curve, i.e., 
JC = 1. Under these conditions, ordinary perturbation 
theory 00 can be applied. By assuming continuity 
of the Bardeen potential and the well known conserved 
quantity £ (defined below), it was then found |l5| that 
the scale invariant spectrum does not survive the bounce, 
with the actual resulting spectrum being much lower than 
the observed one. The temporary conclusion of this fact 
is that in order that the ekpyrotic model be still com¬ 
patible with the observation, a new procedure must be 
applied to the bounce. 

Arguing against GR during the bounce epoch sounds 
natural, as in particular either the real theory is at least 
5-dimensional, or, worst indeed, in the case of the new 
scenario [j9|, the manifold becomes (curvature) singular 
there, obviously leading to a breakdown of ordinary GR 
across the bounce. In this case, a new criterion should be 
derived to replace the ordinary junction conditions. Such 
a criterion was provided in Ref. |kJ , although without a 
physical motivation, leading to the recovery of the ob- 
servationally correct spectrum. The very exhibition of 
junction conditions leading to a scale invariant spectrum 
could then be seen as a hint that constructing a real¬ 
istic theory satisfying observational constraints was not 
impossible. 

Even if one is prepared to accept such drastic changes 
in the standard cosmological picture, one might wonder 
as to the use of perturbation theory on top of an oth¬ 
erwise singular background |lj|. Moreover, it should be 
mentioned that although the old scenario, because de- 
scribable as an effective bounce occurring at a low enough 
temperature, was avoiding the over-production of grand 


unified scales monopoles [Q, the new model, being sin¬ 
gular, poses this problem in a way which is as acute as 
it was before the advent of inflation. Finally, t he p uzzle 
of trans-Planckian scales J|(J, quoted in Ref. J|l) as a 
caveat for inflation, can be transposed in the new ekpy¬ 
rotic model in the same words. 

This article is organized as follows. After a brief re¬ 
minder of the ekpyrotic model of the universe (Sec. II), 
we examine in detail the junction conditions suggested in 
Ref. ]K| (Sec. |I|). We concentrate in particular on the 
fact that this proposed criterion rests on an altogether 
arbitrary (hence unphysical) normalization function, so 
that whatever spectrum can be obtained: obtaining a 
scale invariant spectrum in this model thus turns out to 
be equivalent to imposing it from the outset. We also 
demonstrate that the new prescription leads to an incor¬ 
rect prediction in the exact case of a radiation to matter 
domination transition. 

We then consider a second possibility, i.e., we examine 
an effective bounce in a context where the linear pertur¬ 
bation theory is still valid. We therefore considered first, 
in section IV, the simplest case in which not only does GR 
apply, but also in which all the calculations can be per¬ 
formed analytically and consistently (indeed providing 
a nice textbook example for cosmological perturbation 
theory illustration), namely that of a K, = 1 bouncing 
universe with hydrodynamic perturbations |2^ ]. Then, 
using this toy model, we examine how the relevant per¬ 
turbed quantities behave through the bounce. We pay 
special attention to the “short time bounce limit” (this 
is related to the question “how sharp is sharp” evoked 
in Ref. [[L4j) and study whether, in this limit, the bounce 
can be considered as a surface where the equation of state 
jumps. If so this would allow us to use the standard junc¬ 
tion conditions. 

The second example that one can treat completely is 
by considering a test scalar field. Indeed, in this case, 
one does not need to specify what the origin of the back¬ 
ground evolution is. In section [v|, we calculate the spec¬ 
trum of a spectator scalar field in such a bouncing back¬ 
ground. Assuming no strong deviation from GR at the 
perturbed level (we remind that such deviations are nec¬ 
essary in the bounce region), and K, = 0, this also gives 
the tensor perturbation spectrum. The description of a 
bouncing universe with /C = 0 requires special care, as 
GR does not allow for such a configuration to take place 
unless the Null Energy Condition (NEC) is violated. Al¬ 
though this case is clearly contrived, it provides at least 
an example where some arguments presented recently in 
the literature can be implemented concretely, at the level 
of equations. 


II. THE EKPYROTIC SCENARIOS 


The ekpyrotic model is supposed ||, to stem from 
the theory by Horava and Witten (?]] and some particu¬ 
lar construction of heterotic M—theory j(j]. It finds its 
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inspiration in the extra dimensional scenarios, a la Ran¬ 
dall - Sundrum [|], and can be motivated by compact- 
ifying the action of 11 dimensional supergravity on an 
S 1 /Z 2 orbifold, compactified on a Calabi-Yau three-fold. 
This results in an effectively five dimensional action read¬ 
ing 
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(1) 


where <j> is the scalar modulus, and T the field strength of 
a four-form gauge field. Two four-dimensional boundary 
branes (orbifold fixed planes), one of which to be later 
identified with our universe, are separated by a finite gap. 
Both are BPS states |li|, i.e., they can be described at 
low energy by an effective N = 1 supersymmetric model, 
so that their curvature vanishes. This is how the flatness 
problem is addressed in the ekpyrotic model. 



FIG. 1: Schematic representation of the old ekpyrotic model 
as a bulk - boundary branes in an effective five dimensional 
theory. Our Universe is to be identified with the visible brane, 
and a bulk brane is spontaneously nucleated near the hidden 
brane, moving towards our universe to produce the Big-Bang 
singularity and primordial perturbations. In the new ekpy¬ 
rotic scenario, the bulk brane is absent and it is the hidden 
brane that collides with the visible one, generating the hot 
Big Bang singularity. 


at rest. Non perturbative effects yield an interaction po¬ 
tential between the visible and the bulk brane. The dis¬ 
tance of the former to the latter can be regarded as a 
scalar field living on the four dimensional visible bound¬ 
ary brane whose effective action is thus that of four di¬ 
mensional GR together with a scalar field ip evolving in 
an exponential potential, namely 





IjdJif-VW 


( 2 ) 


with 


v(<p) 


— Vi exp 





( 3 ) 


where 7 is a constant and k = 8nG = 817 / 171 ^. Apart 
from the sign, the potential is the one that leads to the 
well known power-law inflation model if the value of 7 
lies in a given range ]24| ]. 

The interaction between the two branes results in one 
(bulk or hidden) brane moving towards the other (vis¬ 
ible) boundary until they collide. This impact time is 
then identified with the Big-Bang of standard cosmol¬ 
ogy. Slightly before that time, the exponential potential 
abruptly goes to zero so the boundary brane is led to a 
singular transition at which the kinetic energy of the bulk 
brane is converted into radiation. The result is, from this 
point on, exactly similar to standard big bang cosmology, 
with the difference that the flatness problem is claimed 
to be solved by saying our Universe originated as a BPS 
state (see however |23|1). 



V 


In the “old” scenario Q, the five dimensional bulk is 
also assumed to contain various fields not described here, 
whose excitations can lead to the spontaneous nucleation 
of yet another, much lighter, freely moving, brane. In 
the so-called “new” scenario [;||, and its cyclic exten¬ 
sion Q, it is the hidden boundary brane that is able 
to move in the bulk. In both cases, this extra brane, if 
assumed BPS (as demanded by minimization of the ac¬ 
tion) is flat, parallel to the boundary branes and initially 


FIG. 2: Scale factor in the new ekpyrotic scenario. The 
Universe starts its evolution with a slow contraction phase 
a oc (— r/) 1+l3 with /3 = —0.9 on the figure. The bounce itself 
is explicitly associated with a singularity which is approached 
by the scalar field kinetic term domination phase, and the 
expansion then connects to the standard Big-Bang radiation 
dominated phase. 
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III. COSMOLOGICAL PERTURBATIONS IN 
THE NEW EKPYROTIC MODEL 

A. The background 

As mentioned above, although the physics which de¬ 
scribes the evolution and the collision of the branes is 
very complicated, it is assumed that it can be described 
by means of a simple four-dimensional model. In this 
case, the equations that govern the system are nothing 
but the Einstein equations 



where a prime denotes a derivative with respect to the 
conformal time rj. The Hubble parameter can be ex¬ 
pressed as H = H/a where H = a'/a. The equation of 
state u) can always be written as 


B. The scalar held era 

Let us start with the scalar field era. It is assumed 
that the evolution of the four-dimensional background is 
governed by the scalar field potential of Eq. (||). This 
is a well-studied case and the resolution of the Einstein 
equations leads to a solution where the scale factor is a 
power-law of the conformal time 

a W) = to(-v) 1+0 , 0 ) 

<p(v) = Vi + ■rry ^ 1 + P) ln (~v)- ( 10 ) 

As already mentioned, the function 7 is constant and its 
value reads 7 = (2+/?)/(l+/3). In the ekpyrotic scenario, 
one has 0</5+l<l. Since the Hubble parameter is 
given by aH = (1 + /3 )/77 < 0, this corresponds to a very 
slow contraction. 

In the case where a single scalar field dominates, the 
evolution of density perturbations can be described by 
means of a single equation, 




-1 


- 1 , 


( 6 ) 


fi + 


K/tT 


V = 0 , 


(ii) 


where the function T(j]) is defined by 


r = 


K_ K_ 

n 2 + n 2 ' 


where k is the comoving dimensionless wavenumber and 
the quantity f .1 is related to the Bardeen potential Fourier 
(7) component <I> by the following relationship 


This last function is a direct generalization of the quan¬ 
tity 7 = 1 —7f'/7f 2 for non spatially flat universes since it 
gives zero in the case of de Sitter spacetime. In the case 
of spatially flat sections, the equation of state becomes 
u> = 2y/3 — 1. For a constant equation of state, the func¬ 
tion T or 7 are constant. This is the case for the potential 
(^) and the function 7 gets a constant value, explaining 
why we used the same symbol to denote these a priori 
different objects. For the equation of state u = — 1 (i.e., 
de Sitter space-time), they vanish. Finally, the sound 
velocity can be formally defined as 

cl^p'/p'. ( 8 ) 




( 12 ) 


These quantities are related to the functions v and z 
in Eqs. (19) and (21) of Ref. JlfJ through v oc p and 
^ oc a,yj 7 . The initial condition for the function p, are 
fixed by the assumption that the quantum fluctuations 
are initially placed in the vacuum state. This amounts 
to 


lim fi = -■=— 

k/(aH)^+oo m pl v2 k 


(13) 


As already mentioned, in the ekpyrotic universe, it is as¬ 
sumed that /C = 0. As explained above, the pre-impact 
phase consists in a scalar field dominated era and an hy- 
drodynamical era. We will follow in details the evolution 
of the perturbations during these two eras. 


For a power-law scale factor, the equation of 
motion for the quantity /i can be solved ex¬ 
actly in terms of Bessel functions as /i = 
(krjj 1 / 2 [Ai{k)Jp + i/ 2 {krt) + A 2 {k)J_^ +1 / 2 ){kri)\. In 
this case, the Bardeen potential can be written as 


J 


®(v) 


nv 7 

2 ka 


(M 1/2 


Ai(k)J/3 + 3/ 2 (kri) - A 2 (fc)J_( /3+ 3/2)(fcr?) , 


(14) 


where the coefficients Ai(k) and A 2 (k) are given by 


Ai(k) 

A 2 (k) 


nV 8 Q i ( kr K-’ K f 3 / 2 ) 
m pl cos 07 t V 2 k 
-iA^P. 


(15) 

(16) 


Although these formulas are exact, there are not so 
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easy to work with and their interpretation is not espe¬ 
cially illuminating. In order to facilitate the interpreta¬ 
tion, it is interesting to proceed as follows. In the general 
case, i.e., even if the scale factor does not behave as a 
power-law of the conformal time, the quantity /r can be 
expressed as 


where the coefficients b 2n (if) can be found by plugging 
the previous equation in the equation of motion for /r 
and by identifying the corresponding order in k. One 
finds 


Kv) = 62 n{jl)k 2n , (17) 

71—0 


J 


P(v) = A i 


1 - k l 


dr 


(ci 2 7)(t) 


+^.2«V7 


dr 


(« 2 7 )(r) 


1 -k 2 


dr'(a 2 7 )(r') + 0 (k A ) 


dr'(a 2 7 )(r') 


dr" 


(a 2 7 )(r") 


<0(k A ) 


(18) 


Inserting this expansion into the expression of $, one 
obtains at leading order 


®(v) 


A 2 (k ) H 
2 k 2 ^ 


Ai(k) H 
2 ^ 



(19) 


It is necessary to push the expansion to second order 
because we see that the first term of [fj,/(a-y/ 7 )]' vanishes 
in the expression of the Bardeen potential. Notice that 
Eq. (22) of Ref. [ [UJ is not correct since one integration 
is missing, see Eq. (|l£|). The fact that it is necessary to 
push the expansion up to second order to obtain the first 
non-vanishing term in the Bardeen potential (the same 
is true for the quantity see below) has been called a 
subtlety in Ref. [[To| whereas this fact has been known 
for a long time in the literature, see Ref. [|r| ||). One 
can easily check that taking the limit kr] —> 0 in Eq. (0) 
leads to the same dependence in conformal time as in the 
previous equation. Explicitly, one has 


®(v) 


A 2 (k)l + f3 \ — 3 — 2 /j 
2 fc 2 K V) 


Ai (fc) 2 + (3 
2 3 + 2/3' 


( 20 ) 


In an inflationary universe where (3 — —2, the ^ 41 - 
constant mode is the dominant mode since the ^ 2 -mode 
decays as — 77 when 77 —■> 0. The fact that it is propor¬ 
tional to (3 + 2 and therefore vanishes when [3 = —2 is 
just the well-known fact that there is no density pertur¬ 
bations at all in a pure de Sitter phase. In the ekpyrotic 
case, the situation is exactly the opposite, i.e., the ^ 41 - 
constant mode is no longer the dominant mode. The 
dominant mode is now the ^ 2 -mode since this one scales 
like — 1 / 77 . The important point is the ^-dependence of 
this mode. This can be found by comparing the exact 
equation ([lj) with Eq. (|l9|) which allows us to make the 
link between the constants A±, A 2 and A±, A 2 . One ob¬ 
tains 


~ /j-/3-5/ 2 , Ai ( k ) „ fc /3+i/ 2 . (21) 


Therefore, for (3 ~ —1, the dominant mode acquires a 
scale invariant spectrum since oc k ~ 3 / 2 . Note however 
in that respect that the ekpyrotic and de Sitter cases are 
already at this stage very different since in the de Sitter 
case, the scale invariant part of the spectrum is time 
independent, contrary to what happens in the ekpyrotic 
situation. 


C. The hydrodynamical era 

As the bulk brane is approaching the visible brane, 
the scalar field blows up. In the ekpyrotic scenario, it is 
assumed that the shape of the potential changes and goes 
to zero as the collision is taking place. Consequently, just 
before the collision, the equation of state tends toward 
a stiff equation of state, i.e., ojq = 1. More generally, 
Ref. [|loj considers the situation where 

P 2 

iO = — = OJq + U1\T] + IjJ 2 TI + ■ ■ ■ . ( 22 ) 

P 

Being given a general equation of state 01 ( 77 ), it is easy to 
find the corresponding scale factor. It reads 

a{v) =4 exp drjc + ^ j dr'[l + 3 w(t')]| 

(23) 

where £0 and C are arbitrary constants. In this article, 
without any essential loss of generality, we restrict our 
considerations to the truncated equation of state uj = 
1 + 0777 , i.e., ujq = 1. In this case the scale factor can 
be integrated explicitly using Eq. (2.154) of Ref. Q and 
the result reads 


a{rj) = £ 0 


_ V _ 

2(1 + St^o) + 30777 


2/(1+3ojo) 


(24) 
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FIG. 3: Characteristic functions of the new ekpyrotic model 
[See Eqs. (p2|), (pE|), ( pif ) and (pi])]. The parameters are cho¬ 
sen as lo = 1 and = — 1 (the ekpyrotic scenario requires 
loi to be negative). 


To our knowledge, this is a new solution. Moreover, as 
we will show, all the relevant perturbed quantities can 
be calculated exactly in this model. Therefore, this also 
constitutes a new exactly integrable case for the theory 
of cosmological perturbations. In the previous formula 
the constant C has been chosen such that a(0) = 0 in 
agreement with the new ekpyrotic scenario. There is a 
divergence when u\r] = —8/3. This just signals that 
more terms should be included in the Taylor series of 


the equation of state which anyway, under the form of 
Eq. (22), is only valid for in the vicinity of the bounce. 
Moreover, in the collapsing phase, one has uj\ < 0 with 
?7 < 0 so the model is divergence free in the pre impact 
era. It is interesting to see how the physical quantities 
that characterize the model behave. We find 


H{rj) 

np(r]) 

cl(v) 


77(8 


3wi?y) ’ 
48 




5 3 

e + r '” + 6 


3coiri ’ 


(25) 

(26) 
(27) 


where the “sound velocity” is given by Eq. (|j). As ex¬ 
pected, the Hubble parameter and the energy density 
blow up at 77 = 0. These functions are shown in Fig. |] 
On the contrary, the equation of state and the sound ve¬ 
locity are regular. The fact that the sound velocity is 1 
at 77 = 0 can easily be understood. This is a consequence 
of the equation 


w' = -M((*-u>)(l+u;). (28) 

When 77 —0 and u> ^ —1, it is necessary that c 2 behaves 
as c 2 —> u> in order to obtain a finite u/ as TL diverges in 
this limit. 

In the phase dominated by an (effective) hydrodynam- 
ical fluid, the equation that governs the evolution of den¬ 
sity perturbations reads 


<F" + 3(1 + cl)H& + [2 H' + (1 + 3c 2 )(H : 2 - £)]$ + c 2 s {k 2 - /C)<F = 0, (29) 

where we have assumed that there is no entropy production. This equation can also be put under the form of an 
equation of motion for a parametric oscillator. Let us define p and 9 by 


6 ri a\p + pj \ npa z 



(30) 


where T is defined through Eq. (]t|). Of course, in the ekpyrotic case, these equations should be used with K, = 0, 
although we wrote them here in their full generality since we will use them in the regular K. = 1 case in the following 
section. Then, the equation of motion of p can be written under the form of a parametric oscillator equation of motion 




p = 0 . 


(31) 


As for the scalar field case, one can solve this equation perturbatively. The solution now reads 


77 ( 77 ) = Bi(k)9 1 - 
+-B 2 ( k)9 


fc2 f 0^7) f ^'(c^Xt') + 0(k 4 ) 

-k 2 J T dr\c 2 s 9 2 )(r') J 


rv dr 
9 2 


6 2 (t") 


0 (k A ) 


(32) 


In the long-wavelength approximation, the solution is 

* = (k)^ + -B 2 (k)^ r%. (33) 


As expected, at this order, the solution does not depend 
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on the sound velocity. In the above equation, the integral 


can be easily performed. The Bardeen potential reads 


J 


. 6sBi(k) B 2 (k) _ , n . .. 

= ^2 + g ^ ir] )2 [~ 12 ^ 17 ? + %^ ir l) +321n(8 + 3wi77)], 


(34) 


where s is the sign of the conformal time 77 . The Bardeen 
potential blows up as 77 is approaching zero and the linear 
theory becomes meaningless. In Ref. 0, it is argued 
that $( 77 ) should not be used. Instead, it is proposed to 
use the density contrast [pT] e m = Sp/p which is linked 
to the Bardeen potential by the relation 


Then the superhorizon solution of e m can be expressed 
as 


fc 2 £i(/c) k 2 B 2 (k) r At 

2 Ha 2 + 21 -id 1 J l 2 ' 


e m 


fc 2 $ 

pa 2 


(35) 


Explicitly, the solution can be written as 


J 


8 k 2 sBi(k) 


^0 


3 9 

I + -W177+—(wi?7 ) 2 


4 k 2 B 2 {k) 
27 uf 


[— 12loiT] + 9(o;ir ;) 2 + 32 ln (8 + 3 c^i? 7 )] 


3 9 

1 + -W177+ — (wi? 7) 2 


(37) 


in which the limit u>± —» 0, being singular, is not applicable. We see that the variable e m is regular at 77 = 0 because the 
divergence has been canceled by the factor 1 /(pa 2 ). Note however that for a constant equation of state u) = loo, this 
variable is not regular if —1/3 < c^o < 1 - The regularity of e m thus depends on the matter content as the singularity 
is approached. If we expand the above equation around 77 = 0, we find 


e„ = 


8 k 2 sB 1 128fc 2 -B 2 


*8 


9 to? 


In 2 


1 + -wi?7 + 0 (t7 3 ) 


9k 2 sBiU}( 


8£ 2 


+ k 2 B 2 (1 + 21 n 2 ) 


77 2 +0(7 ? 3 ) 


(38) 


This equation is in agreement with Eq. (42) of Ref. @ 
with lu 2 = 0. For k ^ 0, the first 0(rj 3 ) is replaced with 
0(k 2 i) 2 \n\ri\). The Bardeen potential and the density 
contrast are plotted in Fig. |4| In Ref. 0 , the solution 
e m has not been expanded in the basis of the growing and 
decaying mode but has been written as 

e m = eo(k)D(k,rj) + e 2 (k)E(k,ri), (39) 


where D = 1 + 30 ^ 177 /4 + 0(k 2 r] 2 In \r/\) and E = rj 2 + 
0(r/ 3 ). The link between the coefficients of the growing 
and decaying modes B 1 , B 2 and the coefficients eo and 
e 2 of the ( D , E) basis is obvious 


eo(&) 


e 2 (fc) 


8 k 2 sBi 

9k 2 sB\uj 2 

K 


128k 2 B 2 

+ k 2 B 2 (1 + 21 n 2 ) . 


(40) 

(41) 


At leading order in k, the previous equations imply e 2 = 
3egw^ / 8, in agreement with the equation in the last line 
of the paragraph below Eq. (45) of Ref. 0 . being given 


that in the present context the variable w^ of Ref. 0 is 
simply 3ta 2 /8; this means that eo and e 2 are of the same 
order in k. The inverse transformation reads 


•*'(*> = ^( 1 + 21 " 2 )e »-W Q ' (42) 
B2{k) = (43) 


Let us notice that if we want to obtain the other terms of 
the expansion, we need to use the higher order terms in 
the expression (|||) of 77 ( 77 ). For example the first (respec¬ 
tively second) branch next-to-leading order term can be 
expressed in terms of elementary functions and of dilog¬ 
arithm Li 2 ( 77 ) [resp. trilogarithnr Li 3 ( 77 )] functions [[27] . 
An expansion of these functions around 77 = 0 reproduces 
Eq. (43) of Ref. fu}. 

Finally, let us end this section by a discussion on the 
quantity C( 7 ?)- This one is defined by the following equa- 































D. Matching conditions 



V 


FIG. 4: The divergence in the Bardeen potential and the 
gauge invariant energy density perturbation near the bounce 
in the new ekpyrotic scenario. Parameters are chosen as In = 
1 and wi = —1. The Bardeen potential 4> is from Eq. (pj) 
with sBi = B 2 = 1, while k 2 sB\ = k 2 B 2 = 1 for e m in 
Eq. @). 
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(44) 

.(45) 


Let us briefly recall under which conditions, the quan¬ 
tity £ can be considered as a constant. First of all, K, 
should be equal to zero such that the last term in the 
above equation disappears. Secondly, there should be no 
entropy perturbations. Thirdly, only the growing mode 
should be considered and we note that it is crucial to 
discard the singular mode |p5f . In order to obtain the 
explicit expression of </, one can insert the formula j3^ ) 
giving fi in Eq. (^). As is well known [|||, it is neces¬ 
sary to push the expansion to second order because we 
see that the first term of (p/0)' vanishes. One obtains 

C = —Bik 2 I* dTC 2 s (T)6 2 {T) + B 2 + (D(k 2 ). (46) 


For the exact model studied here the integral in the above 
equation can be performed exactly. The result reads 


c = - 


PsB i 

"W L 


27uJir] - 


2 + UiT] 


■ 48 ln(— rj) 


Ik) 


(47) 


Therefore, this quantity has a logarithmic divergence as 
the point where the scale factor vanishes is approached. 
This is in agreement with the analysis of Ref. |uj . The 
divergence is again a signal that the linear theory looses 
any meaning. 


We have at our disposal the solutions in each era. The 
goal is now to join them. The first step is to pass from 
the scalar field era to the hydrodynamical era. Since the 
equation of state can be made continuous at this point, 
one has [a] = [a'] = [a"\ = 0. In this case, the usual 
joining conditions can be used and the Bardeen poten¬ 
tial and its derivative are continuous. This means that 
the growing mode in the hydrodynamical era acquires a 
scale invariant spectrum. In other words, B\{k) ~ k~ 3 / 2 
and B 2 ~ fc -1 / 2 because $ has the same shape in both 
eras. The same applies to the other transition, in the 
expanding regime, from domination by the scalar field 
kinetic term to the radiation epoch. 

Clearly, the non trivial step is how to propagate 
the perturbations through the singularity. We have to 
connect the solution in the pre-impact hydrodynamical 
phase with equation of state ui = 1 + ojfr] with the solu¬ 
tion in the post-impact phase with u> = 1 + ujf 77 , being 
given that 7 ^ uf . A priori, this seems simply impossi¬ 
ble because the theory (a fortiori the linear theory) looses 
any meaning (signaled by the divergence of the scalar cur¬ 
vature and/or of the Bardeen potential): how to perturb 
around a singular background? Even if we are ready 
to accept this, the theory suffers from a serious trans- 
Planckian problem since all the wavelengths become at 
some point smaller than the Planck length J|(J. How¬ 
ever, despite these seemingly insurmountable difficulties, 
Ref. [[To| goes on along the following lines. The fact that 
the quantity e m is regular is used in an essential way. A 
first approach would be to impose [e m ] = [e(J = 0. The 
first condition means eg” = £q whereas the second cannot 
be applied since e^(0) = 3aq/4, which is required to be 
different in the pre- and post impact eras since 7 ^ . 

The fact that the derivative cannot be made continuous 
(contrary to the claims of Ref. (kJ) can be directly traced 
back to the fact that the background is singular. As a 
consequence, one cannot find eg”. Then, a new suggestion 
is given to find the coefficients eg” and eg”. It consists in 
assuming that 

£o 0) = e o ( fc )> 4( k ) = 4{k), (48) 

i.e., one assumes that the energy density perturbation as 
well as the second derivative of (e m — eo D) are continuous 
across the bounce. At this point, we would like to stress 
the following remark: the usual matching conditions stem 
from a well defined geometrical requirement J25|, |2^]. In 
physics, in general, the requirement is made that the 
function and its derivative should be continuous because 
the point considered is almost never considered to be a 
singular point. The situation here is therefore extremely 
special as there does not appear to be any physical rea¬ 
son to enforce any matching conditions, especially on a 
variable which, although finite, is the perturbation of a 
diverging background quantity. Let us however press on 
to assert if they lead, as claimed in Ref. (jTo| , to a unique 
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spectrum: this fact in itself would maybe justify a poste¬ 
riori this choice for the criterion. 

Of course, in the post-impact phase, one is interested in 
the coefficient of the growing mode, i.e., in the spectrum 
of $ perturbations, and not in the coefficients Cq and 


62 ■ This is because the growing mode directly pro vides 
us with the spectrum. Plugging Eqs. (^o|) and ( [HI) into 
Eq. and using the continuity condition (|4^) permits 
to evaluate this spectrum. The result reads 


B>(k) 

B>{k) 


-B<(k) 


1 + 2 1 - 
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<2 


,>2 
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,<2 


,<2 


In 2 
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(49) 

(50) 


Therefore, in the limit of long wavelengths, the dominant 
term is Bf(k) and we have B| > (fc) ~ k ~ 3 / 2 at least as 
long as 7 ^ wf 2 . In this case, the spectrum of the 
Bardeen potential is scale-invariant in the post-impact 
phase. Let us now study the new prescription in greater 
details. It is clear that the function e m can also be written 
as 

e m = eo (k)D(k,ri) + f. 2 ^ f(u 1 ,k)E(k,r]), (51) 

/(wi,fc) 

= e 0 (k)D(k,rj) + e 2 (k)E(k,rj), (52) 

where e 0 (k) = e 0 (k), D(k,rj) = D(k,rj), e 2 {k) = 
t 2 (k)/f{wi,k) and E(k,ri) = f(u>i,k)E(k,r]). Of course 
the choice of the basis has no physical meaning at all 


and we can equally well expand e m (ry) in the basis (D, E ) 
or (D,E). Let us remark that we could choose a more 
general change of basis but in the present context the 
continuity of e m would no longer be guaranteed. In the 
standard case, such a change of basis has obviously no 
consequence on the final spectrum as it should. As we 
are going to show, this is not the case for the new pro¬ 
posal of Ref. 0 L With the new basis, the matching 
conditions of Eq. (|48| ) transforms into 

eo( fc ) = £o (k), e^(k) = e^{k). (53) 

This leads to the following expression for the coefficient 
B>(k) 


B^{k) = -i^B<{k)[ f — 


ser i v "' v/ <_ 


,<2 


,>2 


B<(k) 


TT + 2 ( —— ~ 
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UJ 


i! 

<2 


In 2 


(54) 


Now since the choice of / is completely arbitrary, one can 
always choose f > = wf 2 and f < = w< 2 . In this case the 
first term cancels out and we are left with a spectrum 
B>(k) ~ fc ” 1 / 2 corresponding to a spectral index n s = 
3. More generally, one is free to choose the function / 
as a function of k and in this case the spectrum is also 
completely changed. 


In conclusion, the proposal of Ref. |To]] rests on a non¬ 
physical choice. Choosing the normalization of the func¬ 
tion E such that it leads to a scale invariant spectrum 
seems to be arbitrary. The standard junction conditions 
do not depend on the normalization of the mode func¬ 
tions. So even if one admits that we can somehow pass 
through a singularity, it seems that there is no convincing 
way to find a scale-invariant spectrum. This is probably 
because through a singularity any result can be obtained. 


E. Testing the matching conditions: the radiation 
to matter transition 


In a recent proposal [g9|, it was argued that the junc¬ 
tion conditions advocated in the previous section could 
be expressed in a very similar way to the usual junc¬ 
tion conditions. It is well-known 130 that matching con¬ 
ditions follow from the requirement that [h t j] = 0 and 
[. Kij] = 0 where hij is the metric of the spacelike sections 
and Kij is the associated second fundamental form. The 
question is then: on which surface should these condi¬ 
tions be imposed? The standard answer is to match on a 
surface of constant longitudinal gauge (gauge-invariant) 
energy density denoted e g by Bardeen [[L7| . The proposal 
of Ref. |^9) is to perform the matching on a surface of 
constant comoving (gauge-invariant) energy density de¬ 
noted e m . This is the quantity used above and advocated 
in (10 . 
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The standard junction conditions reduce to the conti¬ 
nuity of the Bardeen potential and (j, [$] = [£] = 0. The 
new ones amount to [<I>] = [7dd* + $'] = 0; see Eqs. (20) 
and (21) of Ref. We have taken the surface layer 
pressure to be zero as it was argued in Ref. Q that 
this does not play a crucial role in the present context. 
However, we will come back to this point shortly. If 
[a] = [a'] = [a"\ = 0 , the two sets of conditions are 
equivalent and both lead to [<E>] = [$'] = 0 as already 
mentioned. On the other hand, there exists a situation 
for which the two sets are not equivalent, namely that 
of a sharp transition, i.e., one for which the equation of 
state u> jumps. In order to discuss the accuracy of the 
new proposal, let us examine the case of the radiation 
to matter transition. The advantage is that the exact 
solution is known and then we can compare whether the 
different set of junction conditions reproduce or not the 
correct result. 


In the radiation to matter domination transition, Ein¬ 
stein equations can be solved exactly and the scale factor 
is given by the following expression |jl^| 



(55) 


For rj < Tfeq, the scale factor is approximatively linear 
in the conformal time and the universe is radiation dom¬ 
inated whereas for 77 3> r] eq it is quadratic in the con¬ 
formal time and the universe is matter dominated. The 
freely adjustable coefficient b = \/2 — 1 is chosen such 
that a(r] = r] eq ) = a eq (note that the different choice 
b = 1 is made in |l^]). The superhorizon solution for the 
Bardeen potential is Eq. (33) which, in the case of the 
radiation-matter transition (|55D, can be written as 



_ br\j 77eq + 1 _ 

b 3 (’n/Veq) 3 (bll/'n e q + 2 ) 3 


R by/l 7eq+ 1 

> ? /, ?eq + 2) 3 




1 

br] / r]e q + 1 . 


(56) 


in which the lower bound of the integral in Eq. (|3^) has 
been chosen to cancel the H/a 2 contribution of the sec¬ 
ond branch. From this expression, it is easy to check 
that, for the growing mode, one has 


$(77 > Veq) _ _9_ 
$(?7 < 77eq) 10' 


easily be obtained and reads 


$<( 77 ) 

$>( 77 ) 


vffjLy + lB}, 

277eqU eq \ 77eq J 3 


48B> 

Veqa 2 q 



-5 


-B 


> 

2 • 


(61) 

(62) 


This result is nothing but the standard result of the in¬ 
flationary cosmology, applied to the radiation to matter 
transition. In the same manner, the quantity £ can be 
calculated exactly. One obtains 

C(77 )=B 2 . (58) 

This result is valid as soon as the decaying mode, not 
taken into account here, had enough time to decay. 

Let us now turn to the piecewise solution for which 
the same situation can also be described by means of the 
following approximation for the scale factor 

a K (v) = (59) 

Veq 

a > (v) = + ) ( 60 ) 

4 V W 

where [a] = [a'] = 0 has been imposed in agreement with 
the background junction conditions. For each region, the 
exact superhorizon solution for the Bardeen potential can 


Similarly, one gets the quantity £ as 

(<(r,) = B<, C(V) = B>, (63) 

where, as emphasized above, the decaying mode is as¬ 

sumed negligible. 

We are now in the position to relate the various quan¬ 
tities of interest before and after the transition. For 
this purpose, let us now apply to set of junction condi¬ 
tions. The standard matching conditions stipulate that 
[<L] = [£] = 0. This amounts to 

B> = B< + ( 64 ) 

B> = B<. (65) 

For a sharp transition having [H] = 0 and [ a"] ^ 0, 
impl ying [p] ^ 0 , the matching conditions proposed in 
Ref. [^9| are equivalent to [$] = [$'] = 0. From the very 
definition ( |44| ) of £, these conditions, together with [w] ^ 
0, implies [CjV 0. Therefore, as announced, the two sets 
of junction conditions are not equivalent. Applying the 

















new matching procedure yields 
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'eq a eq 
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IB, 


( 66 ) 

(67) 


In turn, the coefficients Bf and BJ' are fixed by the ini¬ 
tial conditions at some time 77j. Let us express these coef¬ 
ficients in terms of $; and $[, the Bardeen potential and 
its derivative at some initial time r/i <C r/ eq respectively. 
Before the transition, the result is 


i „ i 


$ < (r ? ) = $ i + -77 i $'--7 7i $( 


( 68 ) 



If one uses the standard junction conditions, the Bardeen 
potential after the transition can be written as 


$>(77) = 


16/ 1 . 10 . 77 ? 

9 1 

-I-$i + -??;$ 

10 1 3 ' 1 


1 + - 
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(69) 


leading to the correct ratio as in Eq. (|57]). On the other 
hand, the matching conditions proposed in Ref. |^9) lead 
to 


$>( 77 ) 
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eq 


1 + -J 

VeqJ 
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(70 > 


The ratio between the constant parts of the Bardeen po¬ 
tential before and after the transition is then very close 
to unity, namely ~ 1 + 77 i<f>(/( 3 $i). Fig. || illustrates this 
point. 

The new proposal does not reproduce the exact result 
for the radiation to matter domination transition. One 
could argue however that the situation for which the new 
junction conditions were suggested is different from the 
standard case, and that therefore new rules must be ap¬ 
plied. This would require different physical prescriptions 
for different situations, whereas it seems to us that a uni¬ 
fied approach is more satisfactory. 


IV. HYDRODYNAMICAL BOUNCE AND THE 
CONSERVED QUANTITY c, 

We now turn to the second topic of this work in relation 
with the bouncing phase. From now on, we shall consider 
a regular bounce, i.e., one such that the scale factor never 
vanishes. As already mentioned, it is clear that such a 
bounce cannot be described in the same framework as 
in the ekpyrotic case because it is impossible to have a 
bounce if AH = 0 within GR. For instance, this can be seen 


FIG. 5: Time evolution of the Bardeen potential during the 
radiation to matter domination transition. The full line shows 
the exact solution (p6[), the dashed line represents the usual 
approximation (B9[) , and the dotted line is obtained using the 
new proposal ([fOl) for the junction conditions. This last ap¬ 
proximation is not in fact constant but, as discussed below 
Eq. ©. only hardly varying at all on the scale shown. In 
order to ensure that initial conditions are identical for the 
three curves, we have set numerically <E>i = <f>( = 1 [see above 
Eq. (|6^)], and used Bi ~ — b 2 a 2 n? (b&i + 8r] e q <E>[)/9?7e q and 
B2 — 3 / 2 ($i + ?7i$[/3) in Eq. ( |5q ) valid in the limit 771 -C r/ e q . 


if matter consists of a single scalar field since Einstein 
equations yield 


P + P=-^(H 2 -H'), (71) 

Kcr 

which shows that p + p should be negative at the bounce 
(TL = 0, H! > 0) even though it is positive definite, be¬ 
ing also given by p' 2 /a 2 . Indeed, it is well known that 
in order to have a bouncing period in a FLRW back¬ 
ground, one must violate energy conditions that classical 
fluids usually do satisfy (as was presented in Ref. JIT and 
whose origin can be in fact traced back to Ref. jHj). If 
one insists on having a flat K. = 0 situation with, say, a 
scalar field alone, one must either use other equations, or 
assume the existence of a singularity. 

The only way to have a bounce in GR with a well be¬ 
haved (NEC preserving) hydrodynamical fluid as the only 
source of energy momentum is in the case of a closed, 
K, = 1 universe. In this case however, as recently dis¬ 
cussed | 2 ^|, one finds that the bounce must be followed 
by an inflationary epoch, thus considerably lowering the 
interest of the model as an alternative to inflation. We 
shall nevertheless study this case as the only fully self- 
consistent possibility. 

In the literature jib|, [TfJ, it was suggested to treat the 
bounce as a GR sharp transition, i.e., to assume that 
the time scale of the bounce is very short and therefore 
that the theory has “no time” to deviate too strongly 
from GR. The continuous and self-consistent GR model 
developed here will allow us to test the validity of these 
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hypothesis, namely that the bounce can be appropriately 
approximated by a sharp transition between the slow con¬ 
traction phase and the radiation era, and that the curva¬ 
ture perturbation £ is continuous in agreement with the 
standard junction conditions. 

Let us now turn to the description of the model con¬ 
sidered in the following sections. We choose the behavior 
of the scale factor around the bounce to be 


a{p) = £ 0 



(72) 


This choice is reasonable since any function describing a 
bouncing scale factor can be approximated by a parabola, 
at least in the vicinity of the bounce. Any other choice 
would thus be equivalent to this one, and the results one 
would get, for instance by including higher order extra 
terms in Eq. (|72j) seen as an expansion, would be quali¬ 
tatively unchanged. Such a behavior for the scale factor 
results from the presence of an hydrodynamical fluid with 
an unusual equation of state. The various physical quan¬ 
tities needed to describe the bounce such as the energy 
density, pressure, equation of state and sound velocity 
are displayed in Fig. |i]. In particular, the sound velocity 
is given by the relation 


c 


2 
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3 


4 + ^0+^0 (-7 
_ \ 7 7o 

r io — 4 + ^2 + 




(73) 


From the figures, one can already see that the bounce 
is rather unlikely to be well described by a sharp transi¬ 
tion which would require a finite jump in both u> and Cg. 
One could however argue at this point that this is due to 
our approximation for the scale factor at the bounce: an 
even scale factor leads to an even equation of state and 
therefore to a transition which cannot be assumed sharp, 
even if the transition duration goes to zero. 

Another, more important, reason to oppose the sharp 
transition treatment of the bounce lies in the following. 
In the neighborhood of the bounce we have 


{P + P)(V = 0) = 4(l-4 
l o \ Vo 


(74) 


which can be easily interpreted. Restoring the usual units 
(with c the velocity of light), the above equation can be 
rewritten as n(p + p) = 2 /£q[IC — ^/(c 2 ^)], 7n which we 
interpret as the curvature scale [= a(rf)/^J\iC\] at the 
bounce, and to = (oPo/c is the physical time taken by 
light to go across the bounce. It makes sense that some 
exotic matter (p + p < 0) is required if the growth of the 
universe is faster than light. Similarly, the sound velocity 
given by Eq. (||) reveals that if 770 < 1, there always is a 
point at which c 2 diverges. This is connected to the vio¬ 
lation of the null energy condition p3] seen in Eq. (74). 



V/Po 
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FIG. 6: Components of the stress energy tensor (top panel) 
deduced from Eqs. (dj) and (Bias functions of the conformal 
time p/po for the scale factor (h3). Shown are the values of the 
background energy density p (full line), the pressure p (dotted 
line) and their ratio u> (dashed line). In this figure as well as 
the following, the parameters have been given the particular 
values = 1 and po = 2. The bottom panel displays the 
sound velocity c 2 as given by Eq. (|73|). 


This has of course important consequences with respect 
to our wish to have a short duration bounce. In particu¬ 
lar, it means that one cannot, in this framework, investi¬ 
gate the short bounce limit for which po —> 0. Therefore, 
we reach the conclusion that the time scale of the bounce 
cannot be made arbitrary short if we want to deal only 
with well-behaved quantities. This provides another ar¬ 
gument against the sharp transition limit. We shall for 
now on restrict ourselves to the case po > 1 , which is 
consistent with our choice of setting K. = 1 in order to 
avoid unnecessary exotic matter. 

Let us now discuss the standard junction conditions. 
For the background, matching by brute force the pre- 
and post-impact phases (i.e., assuming that the bounce 
time scale is negligible) means [H] 7 ^ 0 since H has not 
the same sign before and after the bounce. On the other 
hand, the junction conditions applied to the background 
demand that [H\ = 0. Therefore, it seems that it is 
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already impossible to use the standard GR conditions at 
the background level, as pointed out in This was 
also the reason why, in Ref. Q , a surface layer pressure 
term was introduced such as to allow for a jump in 7 ~i. 

Accordingly, let us admit that we can study the per¬ 


turbative level. It has been shown in Ref. [^5| that the 
well-known cosmological perturbation matching condi¬ 
tions for KL = 0 also holds for 1C = 1. In this last case, 
the quantity (j is not conserved, and it is better to work 
with another quantity £ BST defined by 


Cbst — o 


n 2 


3(1 + u)(H 2 + IC) 

n 2 „ k 2 
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(76) 



FIG. 7: The parametric oscillator potential 8"/9 [See 
Eqs. (|o|) and (|3l])] as functions of the conformal time with 
the same parameter as Fig. [| 


The derivative with respect to conformal time of the pre¬ 
vious quantity reads 


1 C ' 

TT BE 


2 U 2 lA / k \ 
3(1 + lu)(H 2 +IC) 


(H- 1 ®' + $) . 

(77) 


Therefore, we see that, even if K, 7 ^ 0, Cbst is approxi¬ 


mately constant on superhorizon scales and this property 
has been used in Refs, [jl^, [r|. With the exact toy model 
at our disposal, this can be explicitly tested. 


We now turn to the study of the perturbed quanti¬ 
ties. The effective potential 9"/9, see Eq. (|3^), for the 
scalar perturbations is given in Fig. fj]. We assume that 
there is no entropy production and the equations gov¬ 
erning the evolution of Bardeen potential are Eqs. ( |29| ) 
and ([by). Two cases must be studied. In the short wave¬ 
length limit, for which c 2 (k 2 — K.) 0"/0, </ BST is clearly 

not conserved. Therefore, the only case which remains 
to be studied is that of long wavelengths. The latter ap¬ 
proximation can be applied if c 2 {k 2 — 1C) <C 0"/0, where 
k = 1,2, • • •. This is obviously true for k = 1, i.e., a mode 
that cannot be confused with the background. This is less 
clear for higher k modes, and depends on the parameter 
values. In the short bounce limit 770 —> 1 we are interested 
in, one has c 2 (rj = 0) ~ — ( 5 / 3 )(? 7 q — 1 ) _1 — l/ 3 + 0 (r 7 Q — 1) 
and 9"/9{i 7 = 0) ~ —(15/2)(t7o —1) — 1 -F 3 / 2 -FO(t 7 § — 1), so 
that the ratio tends to the fixed value 9/2. In this limit, 
the mode k = 2 also marginally satisfies the long wave¬ 
length requirement. The approximation breaks down for 
k > 3. Since there exists at least one physically meaning- 
full mode for which the approximation is valid, we can 
proceed and use Eq. (Q) for the relevant modes. In the 
case at hand, the integral can be performed exactly and 
the final result can be written as 
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<h(?7, fc) = .Bi(fc) 
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(78) 


As expected one branch is even and the other is odd. One also sees that if 770 = 1 the minimum of the even branch 
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is zero. This makes sense since in this case p + p = 0, 
i.e., the de Sitter equation of state for which it is known 
that <F = 0. The two branches are plotted in Fig. ||. 




FIG. 8: Odd (top) and even (bottom) modes of the Bardeen 
potential [see Eq. (JT^)] through the parabolic bounce as func¬ 
tions of the conformal time with the same parameter as Fig. 

We are now in the position where we can estimate_C BST 
for these long wavelength modes. Expanding Eq. (|7^) to 
leading order, one finds 

(79) 

In the above relation, only the I? 2 -mode of £ appears in 
the first term since the I?i-term is of order k 2 Bi(k), see 
Eq. (lifil) , whereas the corresponding term in the Bardeen 
potential is of order B\(k). Far from the bounce, when 
77/770 1, the first term in Eq. ( ffsl ) tends to zero, while 
the second goes to —B 2/7 on both sides, as can be explic¬ 
itly checked in Figs. p]. Therefore, if we consider a long 
time interval Arj/po 1, this quantity seems to be in¬ 
deed constant. On the other hand, close to the bounce, 
over typical bounce time scales, i.e., Ap/r/o ~ 1, one 
clearly sees in Figs. |ii| that the quantity C BST is not a con¬ 
stant. This is due to the fact that, during the bounce, the 


“growing” and “decaying” modes are of the same order 
of magnitude as shown in Fig. p]. 



v/Vo 



v/vo 


FIG. 9: Two typical configurations for the supposedly con¬ 
served quantity Cbst as functions of the conformal time with 
the same parameter as on the previous figures. Top panel: 
B\ = 10 and £>2 = 1- Bottom panel: £>1 = 1 and B 2 = 10. 


To conclude this section, let us summarize what we 
have learned from the simple toy model used here: (i) 
we have seen that the standard GR junction conditions 
applied to the background are not consistent with a 
bounce since TL has not the same sign before and after 
the bounce , (ii) we have shown that the equation of 
state does not jump at the bounce, (iii) we have found 
that a bounce cannot be made arbitrary short without 
violating the null energy condition p2| | and finally (iv) 
we have noticed that the quantity (j BST is not constant 
on the typical bounce time scale (recall that beyond the 
bounce epoch, our model looses its meaning and should 
be matched to another era). 
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V. A TEST SCALAR FIELD IN THE 
EKPYROTIC UNIVERSE 


the first one (i.e., the scale factor) is even. Then, it is 
easy to show that 


Let us now turn to the last case for which one can ex¬ 
plicitly calculate the various physically meaningful quan¬ 
tities during a regular bounce. We now assume that a 
bounce took place, and consider perturbations of a test 
scalar field in that background. In this case, we do not 
need to specify the origin of the scale factor. As in sec¬ 
tion 0 we question the conservation of what we know 
is conserved in sharp transitions. We can also regard 
the perturbations studied in this section as gravitational 
waves, provided then that K, = 0 and that somehow the 
necessary modification of GR is negligible on the tensor 
part of the perturbations. 

The equation of a test scalar field in a spatially FLRW 
spacetime with a scale factor given by the previous ex¬ 
pression is 




u 



(80) 


As mentioned above, this is also the equation of motion 
of gravitational waves in GR if 1C = 0. The solution of 
this equation possesses two regimes determined by the 
relative contribution of the two terms k 2 and a"/a. The 
transition time ??j(fc) is defined by k 2 = a"/a. In the case 
of a parabolic scale factor, one has 


a 

a 


Vo 


[l + I( 


2 

\Vo) . 


-l 


(81) 


The maximum of the quantity a"/a is I/^q and define the 
only characteristic scale of the problem, i.e., fc max = 1 /»7o- 
Let us define the parameter e by e = fc/fc max , then one 
has 


r S (fc)=±^V2(l-e 2 )- (82) 

The next step is to solve the equation of motion. From 
the above considerations, we see that there are three dif¬ 
ferent regions. In the first region where rj < —rjj(k), we 
only consider positive frequency modes and we have 

wW = ~j= exp [~ik(r] - 77 i)] (83) 

where r/j is an arbitrary initial time. In the second region, 
where —r]j(k) < ij < r]j(k), the solution is given by 

m(v) = B x a{rj) + B 2 a{r]) f -^r-. (84) 

Jo a 2 (T) 

The lower bound of the integral is a priori arbitrary. 
However, it is very convenient to take it equal to zero be¬ 
cause in this case the second branch becomes odd whereas 


(85) 
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FIG. 10: The expectedly conserved quantity p,/a for a specta¬ 
tor scalar field (or gravitational wave production) in a bounc¬ 
ing universe, i.e., Eq. (jS5|), with parameter fixed at £o = 1, 
r/o = 2., B\ — 1 and B 2 = 1. It is clear on this figure that 
this quantity is not actually conserved through the bounce. 



dr 


Vo 


1 0 ft2 ( r ) 240(17) 
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Finally, the solution in the third region where 77 > ?7j(n) 
can be written as 


mi{v) 


c 1 

y/2k 


exp(-ikri) + 


C 2 

V%k 


exp(+iki~i). 


( 86 ) 


Fig. ^ shows the solution (|85|) as a function of the 
conformal time through the bounce. In the usual situ¬ 
ation, the quantity h = p/a is conserved because only 
the growing mode a plays a role. In the present situ¬ 
ation, it is clear that the usually conserved quantity is 
not actually conserved. The reason is that through the 
bounce the odd mode (which is the decaying mode in 
general) a J ' 1 dr /a 2 now plays a crucial role. Therefore if 
we match the bounce epoch to other eras before and after 
the impact, we see that the usual conservation cannot be 
used. 

Since the conservation law cannot be utilized, one has 
to perform the calculation explicitly. Therefore, the goal 
is now to calculate the coefficients C\ and C 2 . Using the 
continuity of the mode function /r and of its derivative, 
we find 
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Ci(fc) 

C 2 (k) 


e zfc(2r/j +77i) r 'l 

-- {-W + ik 9}(-m)[f + if + ik f)(-Vj)W -ikg\(r]j)j, 

e ikVi ( 'l 

- -Tj-jT | - w + i kg] (- ?7j) [/' + *&/] (?7j ) + [/' + ikf] (-r)i)W + ik d\ (Vi ) j, 


(87) 

( 88 ) 


where we have used the short-hand notation / = a(?y) and g = a(?y) dr/a 2 (r). The final result is given by Eqs. (|8i|), 

( |S8| ) where all the functions are explicitly known except the function r/j = rfj(k). Expanding everything in terms of 
the small parameter e one finds 


Cu(fc)= 


4(krj 0 ) £ 


e 2lV2 (4 + iV 2) + 0(e~ 2 ), C 2 (k) = - 


3iir 


2 v / 2 (fc??o ) 3 


+ 0(e~ 2 ). 


(89) 


The spectrum is defined by the following expression 


k 3 P{k ) 


k - 3 K+Vj) 
2tt 2 a(+?7j) 


(90) 


Since the coefficients Ci and C 2 are of the same order in 
k, as shown in Eqs. (^9|) : the power spectrum k 3 P(k) will 
take the form of an overall power-law amplitude times an 
oscillatory function in k. If we parameterize the overall 
amplitude as fc^s -1 for k -C fc max , one has n s = 1, i.e, 
a scale-invariant spectrum. In the other limit k 3> fc max 
where Eqs. (|l]), ( jg8[) cannot be used, the result is obvious 
since the k 2 term always dominates the effective potential 
a"/a in Eq. (|o|): it is n s = 3. We also note that when 
r )o —> 0, the spectrum blows up and this provides another 
argument against a sharp bounce. 

In order to test the dependence of the spectrum in the 
precise form of the bounce, it is interesting to calculate 
the spectrum for another scale factor. We choose 


a(v) =4y 1+ ■ (91) 

This case is treated in Ref. |}2). Let us notice that when g 
is small (i.e., close to the bounce) the previous equation 
reduces to Eq. (|7|) as expected. Then straightforward 
calculations lead to, still for k k max , 

^> = 2 WP 0 ^’ = 

(92) 

In this case, since C\ = — C 2 at leading order, the oscil¬ 
latory part of the power spectrum contributes in a non 
trivial way to the spectral index of the overall amplitude. 
We find n s = 3 which, by comparison with the spectrum 
obtained in Eq. (|8f]), shows that this spectrum is very 
strongly dependent on the actual shape of the bounce. 
This was to be expected since the bounce has already 
been shown not to be a sharp transition. 


The conclusion of this section is that the spectrum 
of gravitational waves (if we accept the trick that, in 
a bouncing universe, the spectrum of a free scalar field 
can be a good approximation of the actual gravitational 
waves spectrum) is in general more complicated than in 
the inflationary case. A first feature is that there exists 
a prefered scale the magnitude of which depends on the 
details of the model. A second property is that, generi- 
cally, the power spectrum acquires superimposed oscilla¬ 
tions due to the fact that, at last horizon entry, the two 
branches contribute equally. Therefore, the shape of the 
spectrum crucially depends on the details of the model. 


VI. CONCLUSIONS 

The conclusions that can be drawn from this work is 
that it seems impossible to apply any known and well 
motivated criterion to pass through a bounce, whether 
regular or singular, in a model independent way as all 
quantities of interest explicitly depend on the details of 
the underlying model. The ekpyrotic model, although 
a potentially interesting alternative to the inflationary 
paradigm, does pass through such a bounce. Therefore, 
if one really wants to calculate the spectrum in the ekpy¬ 
rotic universe then it seems necessary, first, to consider 
a situation where there is no divergence and, second, 
to provide us with the actual (maybe five-dimensional) 
equations of motion during the bounce, knowing that 
these equations cannot be those of GR. 
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